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The well known concept, to reduce the spatio-temporal dynamics beyond instabilities of trivial 
states to amplitude modulated patterns, is reviewed from the point of view of a formal perturbation 
expansion for general dissipative partial differential equations. For codimension one instabilities 
closed analytical formulas for all coefficients of the resulting amplitude equation are given, with 
no further restriction on the basic equations of motion. Both the autonomous and the explicitly 
time-dependent case are discussed. For the latter, the problem of strong resonances is addressed 
separately. The formal character of the expansion allows for an analysis of higher-codimension 
instabilities like the Turing-Hopf instability and for the discussion of principal limits of the amplitude 
approach in the present form. 


I. INTRODUCTION 

Pattern formation in dissipative systems under non-equilibrium conditions is a classical field of physical science 
and has in particular developed from hydrodynamic problems. In addition, this subject has become recently very 
popular in the context of optical, chemical, magnetic, and even biological systems (cf. ref. [Q and references theirin). 
The renewed interest was partially stimulated by developments in nonlinear dynamics. 

Theoretical approaches rely strongly on analytical perturbation expansions, in order to reduce the equations of 
motion to dynamical relevant quantities. To some extent these concepts are limited to a neighbourhood of an instability 
of a simple state. Such concepts have proven to be powerful tools in the context of low-dimensional dynamical systems, 
even with quite mathematical rigour (e.g. [^). For spatially extended systems with a large number of relevant degrees 
of freedom (i.e. the limit of large aspect ratio in the hydrodynamic context) such approaches have been introduced 
based on a multiple scale analysis |Q,Q. They have been applied to a huge number of concrete examples and can 
even be found in textbooks (e.g. [^). In fact the corresponding reduced description beyond simple instabilities, the 
Ginzburg-Landau equation, is very well known, although its properties are extremely complex and are a subject of 
intensive research. 

Like in normal forms for differential equations, the details of the equations of motion, so to say the ’’physics”, is to 
a large extent contained in the coefficients of the reduced equation. Hence it is desirable to have a closed expression at 
hand, to determine these quantities for an arbitrary equation of motion. Such formulas can be partially found in the 
literature but unfortunately certain restrictions on the structure of the equation or the instability are imposed, 

which may limit the applicability of such expressions. It is one goal to fill this gap and to give these formulas for a 
large class of equations of motion including the case of explicitly time dependent systems. Hence for their evaluation, 
which only requires some simple algebraic computations, it is not necessary to perform the multiple scales analysis in 
each concrete case. In this sense our results apply easily to almost all evolution equations discussed in the physical 
context, considerably simplify the explicit computation of an amplitude equation, and show the common algebraic 
structure among the explicit expression for the coefficients as far as they are available in the literature. 

In order to keep the presentation self-contained we review the complete derivation of amplitude equations within 
the well known concept of multiple scales analysis. Although our formulation follows the standard lines, we stress the 
following properties. The expansion can be understood as a formally exact procedure, and no physical assumptions 
(e.g. on scales), sometimes used to simplify calculations, are necessary. In addition, our approach respects the vector 
type structure inherent in the underlying equation of motion, and the final expressions can be understood as scalar 
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quantities with respect to this structure. Hence with a suitable notation the full calculation is by no means more 
involved as for simple model equations. This observation has two important consequences. On the one hand it opens 
the possibility to analyse higher-codimension instabilities starting from physical equations of motion, by going to 
higher orders in the multiple scales analysis. On the other hand one can discuss the principal limits of the multiple 
scales approach in extended systems. In fact we will present one such example, which occurs within the class of 
explicitly time dependent systems. Finally our approach clearly distinguishes between the multiple scales expansion 
and separate approximations for the evaluation of the coefficients, which may be mixed if the concept is applied to a 
definite equation of motion. 

For those readers who are not interested in technical details, section ^ summarises the notation and the results for 
the instability of a single mode, i.e. the explicit formulas for the coefficients of the Ginzburg-Landau equation. A brief 
discussion of the result is supplemented in section IV. The restriction to one spatial dimension is far from being only 
technical, since a formally satisfactory and sufficient general approach for spatially higher-dimensional rotationally 
symmetric systems is still missing, if one disregards the Newell-Whitehead-Segel treatment of nearly one-dimensional 
extended patterns. The technical derivation is reviewed in section g for autonomous systems. As an application 
to degenerated instabilities, the analysis for the Turing-Hopf instability together with the explicit formulas for the 
complete set of coupled amplitude equations is presented in section The peculiarities which arise in the presence 
of explicitly time-dependent equations are analysed in secti on 0 The phenomenon of strong resonances, which 
is associated with this situation, is discussed in section VII. Finally a few remarks on other higher-codimension 
instabilities and the generalisations for the discussion of spatially non-homogeneous situations are given. 


II. BASIC NOTATION AND RESULT 


We consider a physical system being invariant with respect to translations in space as well as in time. In order 
to deal with very general situations we allow for an A^-component real field $(a;,t). A trivial, that means spatially 
homogeneous and time-independent, state should undergo an instability. Without loss of generality we may assume 
that its value is zero, so that the evolution equation reads 


9 $ 




$ e X G R 


( 1 ) 


Here the linear operator C governs the instability and Af denotes all the nonlinear contributions (cf. eqs.(|l^ and 
®)- 


We presuppose that on variation of the system parameters one mode with wavenumber Qc and frequency ujc becomes 
unstable^]. If we measure the deviation from this instability by the small quantity £^, then the linear operator can be 
cast into 


C = + £^£(2) + C)(e4) ^ (2) 

The most general expression for the linear operator at the threshold, which is compatible with the space-time trans¬ 
lation invariance, reads 


. ( 3 ) 

a. 

where the real matrices take the vector character of the field into account^. A similar expression can be also written 
down for the second order contribution, but it is of no special use in the sequel. The (right-)eigenvalue problem of 
the operator (|^) can be solved in terms of Fourier-modes 


^Either the wavenumber or the frequency may vanish. 

2 The formal expression for the linear operator was written down for convenience only, since a lot of evolution equations 
contain derivatives of finite order. Nevertheless the case of an infinite series, i.e. integral operators, is permitted too. It is 
evident from eqs.(^), (|^ and the subsequent derivation, that our approach applies, if the kernel and the moments up to order 
three admit a Fourier transform. 
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( 4 ) 


( 5 ) 

m-.= Y.(ikrko. ■ 

( 6 ) 


OL 


Here the index v numbers the different branches of the eigenvalue problem. The eigenvalues as well as the vector 
part of the eigenfunctions are completely determined by the complex TV-dimensional algebraic equation (^. By 
presupposition all eigenvalues have a negative real part except for one branch v = where the real part vanishes at 

ReA^'''=^(fc) < 0 (|fc| 7^ |gc|), X^''‘\qc) =■ iuJc, =: Uc ■ (7) 

Finally, we need for technical purposes the left-eigenvectors of the matrix 

, (8) 

and especially the eigenvector at the threshold := Vq^/K 

Let us now turn to the nonlinear contributions. Of course they depend on the system parameters, that means on 
e^, too. But to the perturbation expansion only the expression at the threshold £ = 0 and terms of second or third 
order in the field amplitude contribute. Hence if we put 

A/'[^] = A/2 [^] + As ['^ + 

then the most general expressions of second and third order r 

X2[±] = 

a/S 

(c'-n 

^ Im 

AfM = 

oc/Sy 

Imn 

Here the real tensors (0) and (|^) take the vector character of the equation into account^. 

In order to investigate the motion beyond the instability, we expand the solution of eq.(^) in terms of the small 
parameter e, by taking explicitly the amplitude modulation of the marginally stable mode into account 

$(a;,t) =£[u^e*«‘=“+*“‘=*H(ri,r 2 ,...,Ci,C 2 ,...)+M:e-*«‘=^-*“‘^*A*(...)] + + 0{e^) . (14) 

Here the abbreviations 


) 


OisMltf) , 

ead 

I \ 

( ’ dxf^ J 

I'^m 

d°‘ip d'^ip 


dx°‘ ’ dx!^ ’ dx'^ 


UlVmWn ■ 


( 9 ) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Tn := e^t, Cn := 


(15) 


determine the scales of the slowly varying amplitude A. The evolution equ atio n for A can be derived if we require 
that this representation does not contain secular contributions (cf. section HI). This procedure results in the well 
known Ginzburg-Landau equation 


^For gc = 0 a pair of complex conjugate eigenvalues occurs, because of symmetry. Then Uc denotes one of these branches. 
'^In the sequel the symmetry relations 

w} = u} u, w} = ... 


are used. 


3 









(16) 


d d\ I AI2 A 

— Vttit- I A = rjA + rjylpA + D- 


^dT2 " d^2j ' d^f 

Here the convective velocity v and the diffusion coefficient D are given in terms of derivatives of the critical eigenvalue 

dA(‘^‘=)(fc) 


V = Im- 


D = - 


1 


k=qc 


2 dk"^ 

the coefficient of the cubic term is determined by the nonlinearities 

(HclE + A) 


r = 


{V-cWc) 




a afS 

= a/3 

^ ~ a2,J-aa,4 

a/37 

and the coefficient of the linear term is given by the matrix element of the perturbation 

(?c/(27r) dx 


f] = 


iHchc) 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


(24) 


The integral^ just picks out the Fourier component with wave number q^. (.|.) denotes the usual scalar product in 
. For convenience appendix ^ contains the evaluation of these expressions for the Maxwell-Bloch equations as an 
example. 


III. DERIVATION OF THE AMPLITUDE EQUATION 

One inserts the expression (|l^ ) into eq.(|^) and expands order by order taking the notation (I), ®, (I), and 
(0) into account. 

To first order in e only the linear operator contributes. Since the amplitude A acts as a constant at this order one 
obtains 




(25) 


By virtue of the relation L{k)* = L[—k) this expression is nothing else but the eigenvalue equation (||) of the critical 
mode (cf. (|^)). 

To second order in e one gets 


dt 5ti 9ri 


(26) 




[ 1 ] 


* g—iqcX—iuict _ 


4 iy'"+(A 4 [$]) 


[ 2 ] 


^The case qc = 0, in which only the scalar products appear, is captured by this notations as the formal limit qc 0. 
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Here the additional superscript characters indicate at which order in e the expression in the bracket has to be evaluated. 
First let us consider the linear terms. From the identity 

/ \ BA 

(27) 

and the definition (||), the relation 

follows. In general depends on the wavenumber, so that is not an eigenvector of the derivative L'{qc). However, 
as shown in appendix ^ we can consider this special case without imposing any restriction on the validity of our 
results. Although this step is by no means essential, it helps to simplify considerably the following computations. By 
taking the derivative of the eigenvalue equation with respect to k at k = qc and recalling that the real part has a 
maximum at qc, one obtains 


^i<lc)Uc = 

taking the definition ([l7|) into account. Hence eq.(p^) simplifies to 


a$(2) 




[ 0 ] 


u^e 


iqcX-\-iiJc.t 


dA dA 

’ d^i Bti 


^-iqcx-iuct 


dA* dA* 


d^i dri 


(29) 


(A72 [$])'"] . (30) 


The nonlinear term is easily evaluated with the help of eqs. ( 0 ) and ( |l4| ) 

(A/-2[$])P1 = -2^(0)rjA|2 - [g2q,) - 2^w4]^^e2*9<=“+2*“ct^2 _ ^ (31) 

if the abbreviations (H) and are used. We have to check the secular condition (p^, to determine the solution 
Only the modes contribute, because either the wavenumber or the frequency do not vanish. The 

evaluation of this condition yields 


0 = 


d d 


A . 


(32) 


As a consequence the second and third term on the right hand side of eq.(^) vanish. For the solution one obtains 
discarding transients (cf. eq.(p^)) 

= 2rjA|2 + ... ,6, • ■ ■) + (.. .) , 

(33) 

where the constant of integration may of course depend on the slower scales. 

To third order in e the equation of motion (pj) reads 


a$(3) / a$(2) \ 


[ 1 ] 


dt 


dt ) 


fiA r)/4* 

_|_ ^iqcX + lLOct ^-^ I ^-IQcX-lUJct ^-^ 

' a “ 1 “ a 

dT2 dT2 


— £(0)(j)(3) _|_ ^£(3)(j>(2)^ _j_ ^£( 3 ) 7 ^ ^iqcX+iujat _|_ ^^{0)^-iqcX-iuict 

+ |^£(2)^^e*«‘=^+*“'=‘) • A + |^£(2)y*e-*«‘=“-*“‘=‘) ■ A* + (A72[t])'^’ + (A^ait])'^’ 


(34) 


From eq.(^) we have 


, y 


[1] 


no /9R* 

_ ^iQcX+iuJct ^ -iqcX-iuJct^^^_^_ 

— _c^ o ' o 

dri dri 


(35) 


Here and in the remaining part of this section • • ■ indicate Fourier-modes (with wave vector ±2gc, ±3q£or frequency 
±2tUc) iSwc)) which will not contribute to the secular condition. In the same way we obtain using eq.(^) 
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With the identity analogous to eq. 
d' 


dx° 


[ 2 ] 




one has 




d2] 


,dB 

^6 


a-l ^iq^x+iuict 9A a{a 1) 

d^2 2 

'+^ul^t^A {-if , 

dh^~= 


96 

(36) 

1^2 A 

\a—2 iqcX+iuJct'-' ^ 

^ dil 

(37) 

f)2 A 

„tqcX+ruict 

}u,e ^^2 • 

(38) 


For the evaluation of the nonlinear contributions we again keep in mind that only the resonant Fourier-modes have 
to be considered. Then 




[ 0 ]' 


= re*«‘="=+*“‘=‘|y4pA + ^ 

where for the evaluation eq.(|3^) and the abbreviation ( p0| ) have been used. In the same way one obtains 

(A^3[$])[^' = Ae*«‘="^+“‘=‘|Ap24 + + •■• 
using the abbreviation ( P3[) . 

If we now collect eqs.(p^, @), and the evolution equation (R4h reads 


(39) 

(40) 


a$(3) 




d^i dri 


dB* dB* 


d^i dri 


d^2 dT2 J 


dA* dA* 


d^2 dT2 


1 d^X<-''-fk) 

2 


dk"^ 


u^e 


iqcX-\-iuJct 


+ 


k=qc 

(A4[t])‘"’ + {j^smf^ + 


d^A 1 d‘^X^''->{k) 


dk"^ 


y* g — iqcX — iuJct 


k=qc 


d^A* 

'W 




(41) 


where the nonlinear terms are given by eqs.(g^ and (|40|). In addition we have used eq.(|29|) and the analogous relation 

for the second derivative. _ 

Only the terms written explicitly will contribute to the secular condition, for the reasons mentioned above. Eq.(B4) 
yields 


0 = U 


dB dB 

d^i dn 


dA dA 
d^2 dT2 


^2 A 

+ Dj^+vA + r\AfA 

Ot;i 


(42) 


if one takes the abbreviations (H), (H), and ( p^ into account and recalls, that is invariant with respect to 
translations in space and time. 

Finally one has to separate the amplitudes A and B from each other. For that purpose apply {d/dri — vd/d^f) to 
eq.(^). Then 


0 = 


d d 


B 


(43) 


holds, if we use secular condition ( |32D of the preceding order. Since the general solution of this equation is given by 
B = /o(Ci + vTi) + (^1 — UTi)/i(^i + UTi), but B must not contain a secular contribution, one has /i = 0. Hence 


0 = 


d d 


B 


(44) 


and eq.([42[) results in eq.(16). 
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IV. DISCUSSION 


The analysis of the preceding section has shown, that the dynamics beyond an instability of a single mode is 
generically described by a complex Ginzburg-Landau equation, if the underlying dynamics is autonomous. Since the 
solutions of this reduced equation are bounded for ReZD > 0 and Rer < 0 (cf. the motion is correctly described 

at least on time scales t ~ It is not the aim of this article to go into the details of the properties of the amplitude 
equation, which are itself a field of current research. Let us only mention that the change in sign of ReZ? or Rer 
are higher-codimension instabilities, i.e. a kind of Eckhaus instability and the transition from super- to sub-critical 
behaviour respectively. 

Taking symmetries of the underlying equations of motion into account, the apparently complicated expressions for 
the coefficients of the amplitude equation may be simplihed considerably. A frequent constraint in this direction is 
a system, being symmetric with respect to space-inversion. Let us restrict therefore the following discussion to this 
case. As the main consequence of the symmetry the matrix (^) determining the stability is always real. Hence one 
real or two complex conjugated critical eigenvalues for the unique critical mode occurs. 

Consider first the case of a non-vanishing critical wave number ^ 0, which is sometimes called a soft-mode 
instability. Then the restriction to one critical mode implies a vanishing frequency, even for wave-numbers in a 
neighbourhood of Qc- Hence the derivatives of the spectrum and the eigenvectors are real quantities. But then the 
convective velocity vanishes and the remaining coefficients of the amplitude equation (cf. eqs.(^^-(|^) are real, since 
only derivatives of even order occur. 

If the critical wave vector vanishes, the frequency must be finite Uc ^ 0. Such a situation is sometimes called a 
hard-mode instability. Since the eigenvectors are now complex, the coefficients do not reduce to real numbers even 
in the presence of inversion symmetry. But the convective velocity ( 0 ) vanishes for that reason. In addition, since 
Qc = 0, only the spatially homogeneous contributions, i.e. one summand, enter the formulas for the cubic coefficient 
(cf. eqs.(p^-(23)). Of course r reduces to the expression known from the simple Hopf bifurcation (cf. j|, p.l52f]). 


V. TURING-HOPF INSTABILITY 


As stated in the previous section soft- and hard-mode instabilities are typical, i.e. of codimension one, in systems 
which are symmetric with respect to space-inversion. The corresponding codimension two situation, where both 
instabilities occur simultaneously has recently attracted considerable interest even from the experimental point of 
view . For that reason we think it is useful to present the general analytical expressions for the coefficients of the 
corresponding coupled amplitude equations, which govern the dynamics beyond the instability. 

Again we start our analysis from the general autonomous evolution equation (|^) but impose in addition the inversion 
symmetry. Since a soft- and a hard-mode become unstable simultaneously we presuppose that instead of eq.(^ the 
spectrum obeys = 0 and A*'‘^‘=^(0) = iujc- For the weakly nonlinear analysis one has to take both modes into 

account, where the amplitudes may depend on the slower spatial and temporal scales 




(45) 


Since L{qc) is real the vector part uf := of the soft mode is real too. In analogy to section ^ the first order 
leads to the instability condition. For the second order we obtain 


a$(2) 

dt 




- - 


5ri 






-u^e — -e ■= 


5ri 


5ri 


(A4®) 


[ 2 ] 


(46) 


if we again make use of the fact, that the eigenvectors can be assumed to be independent of the wavenumber, and 
that the first derivatives of the eigenvalues vanish by inversion symmetry. That part of eq.(^) which is linear in the 
amplitudes is nothing else but the identical copy of eq.(^) for the two modes. The nonlinear part reads 

(A/'2[$])'^' = -2^(0)rf|A^p -^(2gc)rfe2*«-="^(A^)2 -^(-2q<;)rf*e"2*‘?'=“(A^*)2 

- 2^(0)rf |A ^|2 - [^(0) - 2fw4]rf - [^( 0 ) -b f 

- 2[4(Qc) - - 2[^(gc) + * 

- 2 [ 4 ( qc ) - - 2[Uq^) + . ( 47 ) 

The first two lines coincide with the corresponding expression from the codimension one case (cf. eq.(|M|)). The 
coefficients rf^f, and Ta/b given by eqs.(|^) and (^) evaluated for the soft- and hard-mode respectively. The 
remaining terms describe the interaction of the modes with the coefficient being given by 
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1 


(48) 


r. := - 


^(qc) - iwci ^ 




Since two modes are critical, the secular condition (B4) has to be evaluated twice and yields 

dA^ „ dA^ 


0 = 


0 = 


ari ’ dTi ■ 

The the solution of second order, discarding transients, reads 

^( 2 ) ^ 2rf|24^|2 + +rf + 

+ 2 rf +rfe2*‘^'^*(A^)2 + rf + uf *e“*‘^=*B^* 

_j_ 2P ^iqcX+ii^ct J^S _|_ 2p*g —* A^ * + 2r g“®'3‘=^+®“c*y4^y4'® * -(- 2P*e®'3‘=®““‘=*^^ *j4'^ 

where the constants of integration and depend on the slower scales. At third order we now obtain 

a$(3) 


(49) 


(50) 


dt 




™ S^ax^B^ S-^ax9B^* 
- M* e —-u* e - 


dn 


dn 


S^ax9A^ S-^ax^A^* 


2 dk"^ 

U^q.-^ m 2 dP 


1 • A^+ (£(2)ySg-*9e:r:^ 

H 

tdH * 

^ ^ dn 

e Ac e “ 


dT2 




drn 


k=qc 


q2aS* 

~w 


1 <P\^''‘\k) 

2 


dfc2 




fc=0 


dT2 


dA^ 

— e ° 


2 2 dfc2 




k=0 


9t2 


+ (A/'2[t])'"’ + (A4[t])''’ + • • • 




(51) 


This apparently complicated expression contains the linear part of the two modes which is already known from the 
codimension one analysis (cf. eq.(|^)). The nonlinear contributions read 

(A^lt])'^' = r^e*'?'="^|A®|2A^ + r®*e“*'?'="^|A®|2A^* +r'^e“'=*|A'^P^^ + r'^*e"*“‘=‘|A'f^pA^* 

+ 7^e*'?'=^|A^|2A^ +/*e-*«‘="^|A^|2A^* + 7^e*‘^‘|A®|2A" + 7^*e-“‘|A^I^^^* + ■ •' (52) 

if we restrict the presentation to terms which will contribute to the secular condition. Again the first line describes 
the contribution already known from the analysis presented above, with the coefficients T'^ and T^ being determined 
by eqs.(^^-(^2|) evaluated for the soft- and hard-mode. The new coefficients which mediate the interaction are given 

by 

7^ :=2^(zg,)^C(°^){Af%2rj + 2^(zge)^C(°«{uf,2r:} + 2^(^g,)“C(“°){uf,2rf} {&]R^) (53) 

/3 0 a 

7^ :=2^(-zg,)“(*g,)'5c(“^){uf,2rj + 2^(*g,)“(-*g,)^C(“«{uf,2rj + 2C(°o){uf,2rf} . (54) 

a/3 a^ 

In the same way the other term reads 

_^^S^iq.x\AH^2AS ^gS*^-iq,x\AH^2AS* ^§H^iu.t\AS^2AH ^§H*^-iu:,t^A^^2AH* _^ . . . ^ ( 55 ) 

where the interaction coefficients are 

:=6^(^(7,)^I?(°°^){uf,uf^uf} {€]R^) (56) 
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If one evaluates the secular conditions arising from eqs.(|5l[), (p2D, and 
one finally obtains the coupled set of amplitude equations 


dA^ 

dT2 

dA^ 

Bto 




iS|2 /iS 


H\2 aS 


and separates the amplitudes A and B 

,sdA^ 


dei 


= r,^A^+ r" |A" + s^\A^\^A^ + D 




dA^ 

~W 


with the coupling coefficients being given by 






s^ = 


iv^\l^+S 




(58) 

(59) 

(60) 


and eqs.(^), (|^), (p^, ( |^). The remaining coefficients are of course not changed compared to the codimension one 
case and are given by eqs. ll)-(il) evaluated for the soft- and the hard-mode. 

We have obtained the well known set of coupled real and complex Ginzburg-Landau equations [^, but with 
the general and closed analytical formulas for its coefficients. The simplicity of the whole derivation for this quite 
complicated instability again justifies the slight amount of formalism presented in the preceding sections. 


VI. EXPLICITLY TIME-DEPENDENT EQUATIONS 

Let us now return to the case of a single unstable mode and consider the case that the equations of motion ( 0 ), 
that means the linear operator C and the nonlinear contributions A/”, possess an additional explicit time-dependence 
with period T and frequency Q = 27r/T respectively. Such a dependence may arise in the presence of external time- 
dependent fields. But one should also keep in mind, that a time-dependency may be introduced also by the trivial 
state, if the equation of motion is cast into the form ( 0 ). For the derivation of an amplitude equation the steps of 
the preceding section can be applied. The calculation is quite similar and we use the same symbols, even if they are 
now T-periodic functions of time. However, two main formal differences occur. On the one hand the stability of the 
trivial state is governed by a Floquet- instead of an eigenvalue problem, that means that the eigenvectors itself are 
T-periodic in time and obey 


(t) + (t) = Uk, (t) (61) 

(fc) - = v[!'^*{t)Uk, t) . (62) 

Here L{k,t) is defined by eq.(|^ with being explicitly time-dependent, and the imaginary part of the Floquet- 
exponents is restricted to the first Brillouin zone, Imfi^'^^(k) G [—12/2, Q/2]. On the other hand the matrix coefficients 
in the abbreviations and (^^, which are actually some propagators of the linear equation, are replaced by the 
corresponding quantity of the time-dependent equation. We will present the details below. 

The main physical difference comes from the fact, that one has to pay attention to strong resonances, which are 
known to be crucial already in the corresponding spatially homogeneous situation H- For that reason it is worthwhile 
to consider the explicitly time-dependent case separately. 

We first consider the situation without strong resonances. Let us assume that either the critical wave number Qc 
does not vanish or that the frequency ujc does not obey a strong resonance condition, i.e. 


, Q Q Q 

7 ^ ± —, ± —, ± — 


(63) 


We now proceed along the lines of section HI. It is not necessary to write down all equations again, since almost all 
formal steps are identical. Therefore only those expressions are made explicit which change. We insert eq.(^0) into 
the equation of motion, keeping in mind that the eigenvector is time-dependent, and compare order for order in e. 

To first order the Floquet-equation (H) is reproduced. The equation of motion for the second order is identical to 
eq. (|^) , where we again profit from the fact that fc-independent eigenvectors can be assumed (cf. appendix ^ . Here 
iv of course denotes the derivative of the critical Floquet-exponent with respect to the wave number. For use of latter 
reference we write down the nonlinear contribution explicitly 
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(64) 




al3 




0.(3 


o(3 


where the T-periodic time-dependencies are indicated. Since either the wavenumber does not vanish or the non¬ 
resonance condition (^) holds, evaluation of the secular condition (B5) leads again to eq.(^). Hence the solution 
discarding transients is given by eq.(|3^) (cf. appendix ^ but with and Ff, being replaced by the T-periodic 
quantities 


0/3 


(65) 


1-M. 


=k^0 


ft-T 




a(3 


r,(t) := r dt' 

° 0/3 


( 66 ) 




— 2iu)cT 


ft-T 




0(3 


The indefinite time-integrals have been reduced to definite ones by using the Floquet-decomposition (B7) and the 
periodicity of the remaining factors. 

We proceed to the third order and obtain eq.(^l|) with the derivatives of the eigenvalues being replaced by derivatives 
of the Floquet-exponents. For convenience we repeat the explicit expression for the nonlinear contributions 


(A4[t])'^' =r(t)e*«-=^+*“‘=‘|H|2H + r(t)e"*'^'=^“*“‘^Vp24* + •■• 


(67) 

( 68 ) 

Here r(t) and A(t) are again given by eqs.(^^ and ( |23D evaluated with the time dependent quantities. Transients 
and terms with wavenumber ±2(7c, iS/^c or frequency ±2a;c, ±3tac, which do not contribute to the secular condition 
by v irtue of the non-resonance condition (^^, have been indicated by • ■ •. The evaluation of the secular condition 
(B5) leads to eq.(|4^) and the subsequent considerations are the same as in section HI. Hence we again obtain the 
amplitude equation (p^. But since the integrand in the secular condition is time-dependent, the integral survives 
and the coefficients read 


= Im 


dn^’'<='>{k) 


dk 


k=qc 


1 (i2^('^'^)(fc) 

2 


k—Qc 


V = 


fo (vc(i)lL(t) + A(t)) dt 

fj'(vc(i)luc(i))dt 

9c/( 2^)4/4/^“ (n,(t)e*«=-| 4 ^| 2 i,(t)e*«=-) dxdt 

Jo {Vc{t)\Uc{t)) dt 


(69) 

(70) 

(71) 

(72) 


One might argue that the evaluation of these formulas is in general impossible. But one should notice that one 
only needs the solution of the time-dependent problem ( p^ ) for one period (cf. eqs.(p5|) and (|6^). This is at least 
numerically a very simple task, so that the evaluation even for quite complicated equations of motion can be performed 
on every computer. But also analytical computations are possible, if the Floquet-problem can be handled, e.g. with 
a separate perturbation expansion. 
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VII. STRONG RESONANCES 


From the analysis of the preceding section it is obvious, that in cases of strong resonances additional terms contribute 
to the secular condition. We will discuss in the sequel the implications of each case separately and therefore assume 
gc = 0 throughout this section. 

a. Quartic Hopf-Hard-Mode Instability Consider the case, that the frequency of the marginally stable mode 


n 

Wc — — 


(73) 


Since this condition puts an additional constraint on the instability, this situation can be roughly classified as a 
codimension-two instability. The analysis of the preceding section up to and including the equation of third order 
(cf. eq.(^) is valid. But now the terms with frequency ±3ujc also contribute to the secular condition. Hence for the 
evaluation of the nonlinear contribution (cf. eqs. (|l^, dl^), (p^, (|3|), d^), (|6^) one has to considers these terms too 


(A4 

(A/'s 


+ 7(t) {A*f + 

= A{t)e^‘^-*\AfA + 

+ S{t) {A*f 4 


Here 


(74) 

(75) 

(76) 

(77) 

denote the additional T-periodic coefficients^, whereas • • ■ indicate those terms with frequency ±2ujc which do not 
contribute to the secular condition. Because of the resonance condition ( [T^ the last summands in eqs.([7^ and ( [75| ) 
do not drop in eq.(^). Instead of eq. (® one obtains the amplitude equation, after having separated as usual the 
amplitudes A and B 


7(t) -.= 20 ^ K(t),r:(t)}e-*^‘ 


dT2 ^ d^2 


The additional coefficient is given by 


A = r]A + r\A\‘^A + s {A*y + D 


Jo {vc{t)\lit) + lit)) dt 


d'^A 


fo (llc(i)lMc(i)) dt 


(78) 


(79) 


whereas for the remaining quantities the formulas of the preceding section apply. In contrast to the usual hard-mode 
instability, the amplitude equation does not possess a phase symmetry. Hence the complex phase of rj cannot be 
eliminated, and one has a two-dimensional ’’unfolding-parameter”. 
b. Flip-Hard-Mode Instability Let the frequency obey 


n 

<^c = TT 
2 


(80) 


The corresponding Floquet-exponent is located at the boundary of the Brillouin zone, that means the Floquet- 
multiplier is isolated and takes the value —1. Such a value induces a period doubling bifurcation, which is of course a 
structurally stable bifurcation of codimension one, so that the condition (80) does not imply an additional constraint 


for the bifurcation. If one takes into account that L(0, t) is a real matrix, then the complex part can be eliminated in 
the eigenvalue equations (^) and (Hj) using the abbreviations 


®A pair of complex conjugate eigenvalues occurs = u)c and w^‘^‘=^(0) = —u)c- The branch r'j, leads to the complex 

conjugate expressions. 

^Since qc = 0 only the term a = P = y = 0 contributes in eqs.(p^, (pl|), (|2^), and (p^. 
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u,{t) =: e— Ut) = -Ut + T)€R^ (81) 

vM = -Ut + T)€R^ . (82) 

Here the real vector u^{t) points into the direction of the centre manifold, which in this case is a Mobius strip. As a 
consequence eq.(^ simplifies to 

$(a:,t) = 2 eUc(t)Ar(ri,r 2 ,...,^i,C 2 ,---)H- 

Ar := ReA , (83) 

so that the field is completely determined by real part of the amplitude. This property is a direct consequence of the 
fact, that the centre manifold in the homo gen eous system is a one-dimensional real manifold. It is now quite simple 
to evaluate the consequences of eqs.(|^), (^), an d (|8^ ) for the amplitude equation. To second order one has the 
result ( ^0[ ) and (|^). But to the secular condition (|B5D now both, the second and third summand on the right hand 
side of eq.(^ contribute 


0 = u 


d 


d 


[A + A*) = 2[v 


In addition, if we insert eqs. 


96 9ri 

and (p^ into the definitions 


d 


d 


96 9ri 
and 


Aj^ 

we have 


Lait) = e2-“*r,(t) = / {u^nut')} dt' e R 


jN 


=k=0 


since the evolution matrix is a real quantity. Hence the solution (|^ ) of second order simplifies to 

= rjt) (A + A*)" + u6i) (B + B*) = 4T,{t)Al + 2^(t)B, . 


(84) 


(85) 


( 86 ) 


To third order we already have eq. where all the cubic terms in the amplitude contribute, that means ( [74[ ) and 
( [tsI ) apply. If we evaluate the coefficients (|^, (|^), (|7^), and (ff^ ) using eqs.(^, (|l|), (^), and ( ^ we obtain 

(87) 

( 88 ) 

where 


(Af2[$A])'^’ =f(<)(A + A6^ + -- 
(A/'s [t A])'"’ = A(t)(A + A6' + " 


^(00) 

qooo) 


m = 36"“‘r(t) = 2cr e R 


N 


Ait) = 3e*-“‘A(t) = nr'^’ e IB 


■N 


The secular condition (|B5|) then leads to the amplitude equation 
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92 A^ 


(89) 

(90) 

(91) 


if the amplitude Br is separated as usual. Again the velocity and the real diffusion coefficient are determined by 
the derivatives of the spectrum (|^ ) and (^), whereas the remaining quantities are expressed in terms of the centre 
manifold coordinate 


/o^ (vcit)\tit) + Ait)) dt 

T = --- - - ^ IFl 

Jo iVcit)\uJt)) dt 

<lc/i‘2TT) fo (^it)\A^^\^it)) dxdt 

Tj = ----- 

Jo ivcit)\Mt)) dt 


e R 


(92) 

(93) 


Eq. (0), which is entirely real including the amplitude, should not be mixed up with the real Ginzburg-Landau 
equation. It is sometimes called a Fishers equation 0- 
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c. Cubic Hopf-Hard-Mode Instability Yet there appeared only modifications in the amplitude equations, but the 
general perturbation scheme was not influenced. This feature changes considerably for a third order degeneracy in 
the Floquet-multipliers, i.e. 


Wc = 


n 

3 


(94) 


It is evident from eqs.(30) and (|^ that the nonlinearities contribute to the secular condition. Indeed we obtain^ 


/ d d 


A = a{A*f 


(95) 


with a coefficient a being determined by quadratic nonlinearities. It is not difficult to show, that for almost all initial 
conditions eq.(|^) yields an unbounded, that means a secular, solution. This feature is by no means amazing, since 
a stabilising cubic term, known to be important in the normal form of the spatially homogeneous system, is missing. 
One may cure this defect by tracing back to what is sometimes called a re-summation of the secular conditions Iq, 
p.318f]. Although such an approach seems to be quite common, it is difficult to estimate the validity of this procedure. 
Especially the actual expansion parameter is unclear and one mixes the different scales in an uncontrolled manner. 
In fact, the formal concept of the multiple scale analysis implies, that the secular conditions at each order have to 
be satisfied separately. In this sense the codimension-two bifurcation of this paragraph cannot be treated by the 
perturbation scheme. 


VIII. CONCLUSION 


It was shown by explicit calculation that for every codimension-one instability of a trivial state the dynamics beyond 
the threshold is governed by a Ginzburg-Landau or a Fishers equation. Especially for the frequently met case of 
a soft-mode or a hard-mode instability in autonomous systems, there arises a unique expression for the coefficients 
of the amplitude equation, which covers both the real as well as the complex Ginzburg-Landau equation. The 
reader may object that all these results can be obtained simply from normal forms of ordinary differential equations 
supplemented with symmetry considerations, and that its is not necessary to go through the explicit derivation. 
Alongside the disadvantage, that such an approach is incapable to yield numerical values for the coefficients, e.g. to 
locate transitions from super- to sub-critical instabilities, one has to be careful concerning the validity of such 
phenomenological amplitude equations. It might happen, that in the stage of the derivation secular conditions occur, 
which put severe constraints on the validity of th e mu ltiple scales approach. One such constraint, which in fact 
invalidates the approach, was presented in section VE in conjunction with strong resonances. A more prominent 
example is known in the context of counter-propagating waves |p^ . Here the thorough approach leads to a nonlocal 
reduced description, whereas the validity of the phenomenological amplitude equations containing convection terms 
is limited to a higher-codimension instability. A similar phenomenon is known in the context of phase turbulence, 
described by the Kuramoto-Sivashinski equation Hence symmetry considerations and normal forms of ordinary 
differential equations are very good guidelines for the resulting amplitude equation, but they do not substitute a 
formal derivation. 

The general and formal approach has shown that it is by no means essential to implement certain properties of the 
perturbation expansion with ad hoc assumptions. In fact all scales which seem to be superfluous drop by itself, which 
emphasises the consistency of the perturbation expansion. One should however stress, that the whole physics (or 
mathematics, depending on the readers taste), is contained in the expansion ( p^ ) of the solution. All the remaining 
steps are just straightforward. Hence, by a change of this expansion it is obvious, that different situations, e.g. higher- 
codimension instabilities, phase equations etc. can be also handled for a general equation of motion. Since the whole 
formalism is quite simple, it is indeed possible to treat higher-codimension bifurcations, and to obtain the amplitude 
equation from the basic equation of motion as demonstrated on the example of the Turing-Hopf instability. Different 
cases, e.g. the degenerated soft-mode instability which require perturbation expansions of higher order have been 
already treated ||^ , and will be published elsewhere. 

Finally the formal and general nature of the presented treatment may clearly indicate the principal limits of 
the approach by amplitude equations. Although it was not our purpose to touch the asymptotic properties of the 


®A linear and a diffusive term can be introduced by using a different scaling of the amplitude and the spatial coordinate with 
e. It does not seem to change the subsequent considerations. 
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expansion from the mathematical point of view, our scheme probably contributes to this field as well, e.g. one may 
construct a clear connection to the normal form theory of low dimensional dynamical systems. But a thorough 
investigation will require methods which are beyond the scope of this publication. 


APPENDIX A: TRANSFORMATION PROPERTIES 

Suppose depends on k. We can perform a unitary transformation R{k) depending continuously on k, so that 

= Uc: E{qc) = i (Al) 

holds. The symmetry relation R*{k) = R{—k) may be imposed^. We are now going to use this rotation to transform 
the full partial differential equation (Q) in order to obtain fc-independent critical eigenvectors. The main trick however 
is, that the final results (cf. section y) are scalars so that the transformation cancels in these formulas. Hence the 
evaluation of these expressions based on the original equation is permitted and one may not take any notice from the 
transformation. 

To perform the transformation (^) on each Fourier-mode let us define the real operator 

^ ^ (A2) 


(A3) 


Its inverse can also been introduced since R{k) is unitary. With the transformed field 

$(a:, t) := t) 


the equation of motion reads 


5 $ 

'm 


= ncn-^^ + 7^A/'[7^■^$] =: £$ + ^[$1 


(A4) 


We now apply the derivation of section III, since by definition (Al) the eigenvector of 

Uk) ■■= R{k)L{k)R~\k) 


(A5) 


does not depend on k. Hence we obtain the result of section || but of course the coefficients (0), dl^, and 
are expressed in terms of the new quantities. We have to show that the transformation drops from this expression. 

First the eigenvalues X^'^''{k) are independent of the transformation TZ, so that expression 0 ) and (^ can be 
evaluated in terms of the original quantities (^). Second the transformed field (^) has an amp litude A which may 
differ from the amplitude A of the original field. But if we insert the definition (^^ into eq.(|A3), we clearly observe, 
taking eq. into account, that both amplitudes coincide up to order e. Third the transformation drops from the 
matrix element (p4|), if the definition for C (cf. eqs.(|A^) and ([A^ )) is inserted. Finally we have to show, that the 
nonlinear coefficient can be evaluated from the original quantities. 

For that purpose consider the transformed nonlinearity A^, insert a field consisting of two Fourier-modes '0 = 
0^ exp(z( 7 ix) + 0^ exp{iq 2 x) + c.c., and equate the Fourier-component exp(i(( 7 i + q 2 )x). Then, according to the 
definition 0) and the transformation (|A^) one obtains 


^(i9i)“(ig2)^C^“^^ = 0(%)“(*<72)^^(gi +92)C'(“^^ 

a/3 a/3 


(A6) 


This formula tells us how the tensors change under the transformation. We now repeatedly use this relation to evaluate 
the nonlinear coefficient by choosing 9i, < 72 , 0^ , and 0^ appropriately. Consider the definition of F^ (cf. eq. 0)) and 
apply relation ([A^). Then 


®We need the transformation as a formal tool only in a neighbourhood of |fc| = \qc\. Hence we need not worry about the global 
continuation. 
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r = - 


='' ^ a/3 


holds, where eqs.(Al) and (A5) were used in the last step. In the same way one obtains 

= Ei‘^qc)I.b ■ 


(A7) 


(A8) 


If we now insert both expression into the definition ( P0[) and use again the transformation (M) appropriately we end 
up with 

r = {E-\qc)u„2E-\0)ta}+‘2J2(‘^^qcf{-^qcrmc)C^‘^^^ {g-^-qc)^:,E-\2q,)t,} 

a a/3 

= r . 

In the same way one obtains 


(A9) 


A = A 


(AlO) 


so that the nonlinear coefficient r can be evaluated from the transformed as well as the original equation using the 
definitions dil), (§), (p^, and in terms of the original quantities (0) 

The same approach can be used in the explicitly time-dependent case as well. In that case one chooses a T-periodic 
transformation to make the Floquet-vectors fc-independent 


Eik,t)ul^‘'\t) = u^{t), E{qc,t) = l ■ 


(All) 


The only difference to the preceding considerations results from the time-dependence, which contributes an additional 
term to the transformed linear operator (cf. eq.(A4)) 


Ct — ^ d" 




(A12) 


But it is exactly this property which ensures, that the propagators (^), occurring in the definitions ( |65|) and (|66| ) of 
r^(t) and r{,(t) respectively, transform according to 

g„{t,t')=E{k,m^{t,m{k,t') . (A13) 

Hence the transformation properties (@, (^), (@, and (|A10|) are valid in the time-dependent case too. 


APPENDIX B: SECULAR CONDITION 

Let us determine the solution of the linear equation 

^ = + , (Bl) 

presupposing that the linear operator has only stable modes except for one (cf. section |I^ and that the inhomogeneous 
part, bounded in time, is given by a finite sum of Fourier-modes. 

By considering each Fourier-component separately, the partial differential equation reduces to ordinary differential 
equations and the general solution is easily written down in terms of the matrix (^) 

'(/'(x, f) = ^ f exp [L(fc)t'] Wj,(t — t') + • • • . (B2) 

k ~ 

Here B denotes the constant of integration, and • ■ • that part of the homogeneous solution, which corresponds to 
stable eigenvalues and leads to a transient only. Consider first the non-critical summands, that means |fc| ^ \qc\. 
Then the integrals converge in the long time limit, since all eigenvalues have a negative real part. For the marginally 
stable wavenumber k = qc one eigenvalue with a vanishing real part occurs, so that the integral may increase in time. 
By using e.g. the spectral decomposition of the corresponding matrix 
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(<vr) 


(B3) 


it is evident from eq.(B2) that a secular contribution increasing linearly in time is avoided if the relation 

f-e 

—iujct’ ( 


1 


0= lim ^ 

0—voo 0 




(B4) 


/o 


is fulfilled]^. In this case the solution in the stationary state is obtained by extending the upper limit of the time 
integrals in eq. (|B2| ) to infinity. 

For a T-periodic e xplic itly time-dependent operator similar considerations apply. The secular condition is 
e.g. obtained from eq.(pl|) by multiplication with the critical left-Floquet-eigenfunction (cf. eq.(^)). One 

has to require 


0 = lim — 

0—»oo 0 


^ — iujct' 


{Vc{t')\2AAt'))dt' 


(B5) 


in order to exclude linearly in time increasing contributions. The solution of eq.(Bl) is again determined by the 
expression (B2) with the matrix-exponential being replaced by the corresponding evolution operator — t') of 

the time dependent system. The latter is determined by 


dgAA') 

dt 






(B6) 


The solution in the stationary state is again obtained by extending the time integrals to infinity. This indefinite 
integrals can be reduced to definite ones by taking the Floquet-decomposition 


gAt + T,t')=K,gAt^t') 


(B7) 


into account. Here the constant matrix M determines the Floquet-multipliers. 


APPENDIX C: MAXWELL-BLOCH EQUATIONS 


For purely pedagogical purpose we present the evaluation of the expressions given in section || for the Laser 
instability, which is governed by the Maxwell-Bloch equations [p^ . Thus every formal step is performed explicitly. 
For the physical details however the reader should consult the literature. 

The basics equations of motion, which govern the evolution of the complex valued envelopes e and p of the electric 
and the polarisation field as well as the deviation of the inversion n from the pump level r, read in dimensionless form 


de . 


cr{e-p) 


dp 

dt 


-(1 + zQo)p + (c — n)e 


dn 

dt 


—bn + {e*p + ep*)/2 


(Cl) 


where x denotes the direction perpendicular to the beam. Here the parameters cr, Qqj and b denote the cavity damping, 
the detuning, and the decay rate of the inversion in units of the dephasing rate. The parameter a scales the diffraction 
term. In the sequel we consider Qq < 0 and analyse the instability of the trivial state e = p = n = 0. Introducing the 
five real quantities via 


e = e + z$ 2 ), p = e *^“‘($3 + *$ 4 ), n = $5 


(C2) 


^°Often such a relation is called a solvability condition for the linear eq.(Bl). 
function spaces on which the operators are considered. 


But such a term requires the specification of 
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the equation of motion is cast into the form (|^) 


/ -cr 

1 

0 

1 

a 

0 



( ^ \ 

ilo + adl 

—a 

0 

a 

0 


0 

r 

0 

-1 

0 

0 

$-t- 

-$1$5 

0 

r 

0 

-1 

0 


-$2$5 

\ 0 

0 

0 

0 

-h) 

V $i$3 + ^'2$4 / 


(C3) 


Since the nonlinearity is quadratic, the only nonvanishing contribution using the notation and respecting the 

symmetry constraint reads 

^ (0,0, -uiv^ - u^vi, -U 2 V 5 - U 5 V 2 , 

U 1 V 3 + U 3 V 1 + U 2 V 4 + U 4 V 2 )'^ ■ (C4) 

The matrix is obtained from eq.( p3D by replacing the differential dx by ik. Its characteristic polynomial reads 
0 = [A + 6] [nl + a^(l - rf + 2X{nl + (1 - r)(l + a)a} + A2{(1 + af + ^1+ 2a(l - r)} + 2X^1 + a} + A^] (C5) 


with the abbreviation := Slg ~ To determine the instability threshold r = Vc one inserts A = itUc and obtains, 
by separating the real and imaginary part, the critical frequency 


UJ 


2 

c 


+ (1 - rc)(l + g)g 

1 + cr 


(C6) 


and the threshold 


= 1 


^qc=0 

1 + cr 


= 1 + uj: 


(C7) 


The instability occurs at = 0 since the threshold (C7) is minimal for that wavenumber. The right and left 
eigenvectors (cf. eqs. (i and (^)) corresponding to the eigenvalue iujc := + a) can be read off from the matrix 

as 


'Y' u u 

Me = -*Wc),l - *Wc,0) , = ——,0) 

1 + lOJc 1 + 


(C8) 


Introducing the deviation from the threshold by r = + £^<5, the perturbation has only two nonvanishing 

matrix elements and eq.(E^ yields 


. (Mc)3(Me)l + {v*c)4{Uc)2 ^ . 
(McIMc) 


1 + cr + iWc(l — cr) 


(C9) 


Keeping in mind, that L{k) depends solely on it is obvious that the convective velocity v vanishes. The second 
derivative, i.e. the diffusion constant, is easily expressed in terms of the eigenvectors ( p 8 |) by tracing back to the usual 
Schrodinger perturbation expansion 


n (Mc)i(Mc)2 - (Mc)2(Mc)i 

U = —a - - —j —7 -= la 

(McImc) 


1 + iujc 


1 + cr + iujc{i — cr) 


(CIO) 


Because of eqs.(p^ and (|C^ ) vanishes identically, and w*} has only one nonvanishing com¬ 

ponent 


C'(°°Hmc,m:} = (0,0,0,0,2)^ . 

Then by virtue of eqs. ( 0 ) and (^ ) we have U = 0 and 

" I (0,0,0,0,1)' 


(Cll) 


(C12) 


so that eq.(| 20 |) reads 
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r = 2c(°o)K, 2r j = - (0,0, -(mJ2,0)^ . 

Since A = 0 because of the absence of cubic nonlinearities we obtain finally for the nonlinear coefficient (|l 


(C13) 


4 -(u;)3(Mc)i - fec)4(Mc)2 
b {Vc\2Ac) 


6 1 + cr + iuJc{l — cr) 


(C14) 


The formal simplicity of this example suggests, that the explicit evaluation is not at all a difficult task, even in 
cases which are usually believed to involve lengthy calculations like the Rayleigh-Benard problem. 
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